Conductance distribution in nanometer-sized semiconductor devices due to dopant 

statistics 
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We show that individual dopant atoms dominate the transport characteristics of nanometer sized 
devices, by investigating metal semiconductor diodes down to 15 nm diameter. Room tempera- 
ture measurements reveal a strongly increasing scatter in the device-to-device conductance towards 
smaller device sizes. The low-temperature measurements exhibit pronounced features, caused by 
resonant tunneling through electronic states of individual dopant atoms. We demonstrate by a 
statistical analysis that this behavior can be explained by the presence of randomly distributed 
individual dopant atoms in the space charge region. 

PACS numbers: 85.35.-p, 73.30.+y, 61.72.-y, 71.55.-i, 



INTRODUCTION 



In semiconductor physics, the influence of doping is 
generally accounted for by a homogeneous shift of the 
Fermi-level, caused by the introduction of free carriers 
in the semiconductor. However, when the dimensions 
of a device are small compared to the average distance 
between individual dopants, the discrete nature of dop- 
ing must be taken into account. Each individual ion- 
ized dopant introduces a Coulomb potential well in the 
semiconductor, locally distorting the potential landscape. 
When the number of dopants in the volume of semi- 
conductor that determines the transport characteristics 
of a device gets small, these random potential fluctua- 
tions cause atypical behavior of semiconductor devices 
Q . This is commonly viewed as one of the fundamental 
limits in the ongoing size-reduction of CMOS-technology 
0, El • Mapping the positions of individual dopants jj] 
and the potential fluctuations they induce 5, ;6] has been 
performed experimentally with various techniques. Fur- 
thermore, the influence of statistical fluctuations due to 
random dopants on device behavior has been subject of 
simulations |7j. 



In this work, we experimentally investigate the effects 
of the discreteness of doping on the transport properties 
of small diodes by comparing many identically prepared 
devices. We find that statistical fluctuations caused by 
randomly positioned individual dopant atoms do not av- 
erage out for very small devices. In contrast, fluctua- 
tions dominate the electrical transport properties of the 
smallest devices and cause large differences in the con- 
ductance of nominally equal devices. Furthermore, we 
demonstrate that at low temperature the Coulomb well 
of a single dopant gives rise to a resonant tunneling chan- 
nel. 



MEASUREMENTS 

In order to allow for measuring the transport charac- 
teristics of many identically prepared diodes, we use self- 
assembly methods to fabricate epitaxial CoSi2-diodes. 
The tip of a scanning tunneling microscope (STM) is 
used to characterize and access the devices individually. 

All experiments are performed in an ultra-high vacuum 
(UHV) system with a base pressure of 5 • 10~ n mbar. 
Self assembled CoSi2-islands are grown on Si-substrates 
(resistivity 0.015 f2cm, doping concentration around 2 • 
10 18 cm -3 ) by evaporating a sub-monolayer of cobalt 
onto a clean 7 x 7-reconstructed Si(lll) surface, followed 
by an anneal at 800 °C for about 5 min. The resulting 
hexagon-shaped epitaxial CoSi2-islands have heights in 
the range of 2-10 nm and diameters in the range of 15- 
80 nm. The inset of Fig. [5] displays an STM-image of a 
typical island. Each island is regarded as a nanometer 
sized epitaxial metal-semiconductor diode. The inter- 
island distances are roughly ten times larger than the 
island diameters. To minimize the effect of surface re- 
lated transport channels, the 7x7 surface reconstruction 
surrounding the islands is destroyed by exposing it to 
atomic hydrogen for 10 min at a substrate temperature 
of 400 °C. 

Current-voltage (IV) measurements are performed by 
positioning the STM-tip over an island. After switch- 
ing off the feedback loop, the tip is lowered by 15 A, 
which is found to be sufficient to make contact to the 
island. Then the current is measured while ramping the 
voltage, yielding well-reproducible iV-curves that reflect 
the properties of the metal-semiconductor contact. The 
sample preparation and measurements techniques are de- 
scribed in more detail elsewhere I8J. 



Room temperature measurements 

Some typical room-temperature measurements are dis- 
played in the inset of Fig^ showing weakly rectifying 
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FIG. 1: The zero bias conductance per unit area as a func- 
tion of island area at room temperature. For large islands, 
the measured values fall within a narrow range. For smaller 
islands, the scatter is rapidly increasing and is much larger 
than the measurement-inaccuracy. The inset shows some typ- 
ical ZV-curves. The numbers in the inset correspond to the 
numbered data points in the main figure. 



TV-curves, compatible with the p-type substrate-doping. 
The overall shape of the curves is very similar. To in- 
vestigate the dependence of the diode's conductances on 
their size, in Fig. ^ the zero-bias conductance per unit 
area is plotted as a function of island area for more than 
40 different islands, grown on two similar p-type (boron- 
doped) samples. One would expect that devices of equal 
size yield the same value for this quantity. Indeed, for 
larger devices the conductance per unit area falls within 
a narrow range. However, towards smaller areas there is 
no definite behavior. Instead, the scatter in the measured 
values increases rapidly and nominally equal devices yield 
very different results. 

To study the increasing scatter in more detail, the stan- 
dard deviation of the conductance per unit area is plotted 
in Fig. |2 as a function of diode area, calculated from the 
data of Fig. ^ We clearly observe that the standard de- 
viation increases with decreasing island size. 

Measurement inaccuracies are not the source of this 
random scatter, as proved by the low noise level (typically 
less than 10 pA around zero bias) and good reproducibil- 
ity of measurements on the same device. Moreover, this 
scatter is absent in similar measurements on low-doped 
samples (not shown). Schottky barrier inhomogencitics 
3, which have been observed at non-epitaxial interfaces 
lOL fill H^ . do not play a role here, as they are due 
to variations in the atomic arrangement at the metal- 
semiconductor interface itself. The metal-semiconductor 
interfaces in our devices are perfectly epitaxial and mono- 
crystalline, which is supported by cross-sectional trans- 
mission electron micrographs of similar structures |l3| . 



FIG. 2: Standard deviations of the measured conductance per 
unit area for various values of the device area A. Each point 
is calculated from eight neighboring data points in Fig. Q 
The dashed line is a least square fit of the function C/y/A to 
these points, yielding C = 1.8 nA/V/nm. The arrow indicates 
which island size corresponds to an average of one dopant 
atom per island ({N) = 1) according to our analysis. The 
inset shows an STM-image of a typical island, acquired at 
room-temperature directly after preparation. 

This also rules out interface defects or grain boundaries 
as the origin of the fluctuations. In addition, BEEM 
measurements of CoSi2-films on undoped Si(lll) reveal a 
perfectly homogeneous Schottky barrier height 11] . Only 
islands showing an atomically flat and defect-free surface 
in the STM-images were included in the analysis. 

We will show that the increased scatter in the mea- 
sured data is caused by the presence of randomly dis- 
tributed dopant atoms in the Schottky barrier. An ion- 
ized dopant atom locally distorts the barrier, giving rise 
to a local barrier reduction and thus a high conductance 
spot (see the inset of Fig. |3J . Then, the observed spread 
in Fig.|2]is directly related to the spread in the number of 
dopant atoms N in the barrier of a device. From Poisson 
statistics, the relative spread is given by 

sd(N) _ 1 

where sd(N) is the standard deviation of N. Because (N) 
is proportional to the island area j4, this gives that the 
standard deviation of the number of dopants per unit 
area increases when the device area decreases. Indeed, 
Fig-0shows that a function of the form C/y/A describes 
our observations appropriately. 

Low temperature measurements 

To further investigate the influence of individual 
dopant atoms on the transport properties of small diodes, 
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FIG. 3: Measured 7V^-curves of two different CoSi2-islands 
on the same sample. The upper graph shows three measure- 
ment of one island (island A, A = 1300 nm 2 ) and also three 
measurements of another island (island B, A — 1500 nm 2 ), 
demonstrating the reproducibility of the measurements. Fea- 
tures occurring at \V\ < 0.1 V (|/| < 1 pA) are due to noise. 
The lower panel shows the normalized differential conduc- 
tance of the same two islands, calculated from an averaged 
TV-curve. Arrows indicate the main features in the in- 
curves, showing up as peaks in the lower panel. The inset 
schematically shows the band diagram of an unbiased device. 
The dashed line is the initial conduction band (Ec) profile. 
The solid line is Ec perturbed by the dopant's Coulomb well, 
causing a local barrier lowering Acp. 



we performed similar experiments in a low-temperature 
UHV-STM operating at 4.5 K. In the upper part of 
Fig. |3 several TV-curves of two different islands are plot- 
ted. Each represents a full and independent measure- 
ment cycle (stop scanning, make contact, ramp V, re- 
lease contact, resume scanning), demonstrating their re- 
producibility. The current is plotted on a logarithmic 
scale, to make the features of the curves visible over sev- 
eral orders of magnitude in current. Note that the noise- 
level is below 1 pA. 

From the measured TV-curves, the quantity (V/I) ■ 
(dl/dV) (the so-called normalized differential conduc- 
tance) has been calculated. It was calculated from the av- 
erage of several TV-curves per island and plotted for the 



same two islands in the lower part of Fig. [3] This is a well- 
known approach in the field of scanning tunnelling spec- 
troscopy and is very useful as it reduces the overall 
exponential behavior and enhances bias-dependent fea- 
tures. Indeed, the weak features in the TV-curves turn 
into clear peaks and the peaks show a one-to-one corre- 
spondence to the features in the TV-curves. As can be 
seen from Fig. |3 features appear at one bias polarity, 
only. This is true for all acquired curves on the same 
sample. For substrates with the opposite doping type, 
the peaks occur at the opposite bias polarity (not shown) . 

The peaks can be explained by resonant tunneling 
through a discrete energy level of a dopant atom, oc- 
curring when the Fermi-level at either side of the barrier 
lines up with an energy level of the dopant's potential 
well. The resulting resonant channel produces a feature 
in the measurement curve. Resonances are expected at 
bias voltages from zero upto the barrier height (roughly 
0.5-0.7 V, depending on the dopant type). This is con- 
sistent with the observations. The actual bias voltage at 
which a resonance occurs is predominantly determined 
by the distance of the dopant atom to the interface. The 
number of peaks in a typical spectrum (1 to 4) corre- 
sponds to the expected number of dopants in the devices, 
as we will show later. This confirms our hypothesis that 
individual dopants influence the conduction path. 



ANALYSIS 

Motivated by the foregoing observations, we present 
a model that links the conductance fluctuations ob- 
served at room temperature to (random) dopant posi- 
tions. Moreover, the parameters in the model are directly 
related to device parameters. 

Because the dopants in the substrate are randomly 
distributed, the number of dopants N in the barrier of 
a given island is Poisson-distributed with parameter A. 
This means that the probability P\(k) that N equals k 
for a certain device is given by 

Px{k) = ^e-\ (1) 

where A = (N) is the mean value of the number of 
dopants in the barrier of this particular device. The pa- 
rameter can be expressed as A = AtN^, where A is the 
area of the island, t the effective thickness of the barrier 
and the average doping concentration in the barrier. 
As shown in Fig. the scatter in the number of dopants 
in the barrier per unit area satisfies sd(N)/A cx \j\f~A 
and therefore is large for small islands. 

We will now show that this effect also creates scatter in 
conductance measurements. We assume that dopants lo- 
cated in a certain region close to the metal-semiconductor 
interface induce a local barrier lowering that gives rise 
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to a low-resistance transport channel 17]. The conduc- 
tance of a transport channel induced by a dopant atom 
depends on the dopant's distance to the interface. Since 
the distance is a random variable, the conductance G\ 
of a single channel is not the same for all channels, but 
has a certain probability distribution. The distribution 
is given by a probability density function fi(g), meaning 
that fi(g)dg is the probability for a certain channel to 
have a conductance Gi between g and g + dg. Both the 
distribution of the position of a dopant in the barrier and 
the dependence of the conductance of a channel on that 
position are contained in fi(g). 

Assuming that the values of the conductance of the in- 
dividual channels are independent and characterized by 
the same distribution given by fi(g), the total conduc- 
tance of a device is given by the sum of the contributions 
of the individual channels. Here, we neglect the back- 
ground conduction and assume that the conductance is 
dominated by these channels. If there are k channels, the 

total conductance is Gk — J2i=i @i \ where the conduc- 
ts) 

tance of the individual channels is denoted by G\ . The 
density function fk{g) of Gk can be calculated explicitly 
by taking the fc-fold convolution of fi(g) 

fk(g) = h{g) * ■■ ■ * h{g) ■ 

" v ' 

k times 

Finally, taking into account the Poisson distribution of 
the number of channels (Eq. the density function 
f(g, A) of the total conductance G of a device can be 
computed as 

oo 

/(<?, A) =£>a (&)/*(<?)■ (2) 

fc=0 

In other words, for given device parameters A and 
t, f(g, X)dg is the probability that the total conductance 
G due to dopant atoms in the diode is between g and 
.9 + dg. 

To illustrate its behavior, Fig.0]shows a map of f(g, A) 
(from Eq. [2J as a function of g/X (proportional to the 
conductance per unit area) and A (which is proportional 
to A). For a fixed value of the device area (that is a 
fixed value of A and a vertical line in the plot) the color 
scale gives the probability density to find a device with 
a particular conductance per unit area g/X. For large 
values of A, all density is concentrated around g/X = i 
while for smaller A it is spread over an increasingly wide 
range of values. In this figure fi(g) was chosen as a 
uniform distribution. However, as we will show next, the 
general properties of /(<?, A) are not strongly dependent 
on the particular choice of fi(g). 

Without making any assumptions on the choice of 
/i(5)i we can compute the moments of f(g,X) in terms 
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FIG. 4: Map of f(g,X) as defined in Eq. ||5J plotted as a 
function of g/\ (vertical axis) and A (horizontal axis). The 
function values are represented in a linear gray-scale, where 
black corresponds to and white to 0.65. For fi(g) a uniform 
distribution was chosen. 

of those of fi(g). First, the mean value of G satisfies 

p OO r. 

(G) = I 5/(5, A) dg = £ P A (AO / gfk(g) dg = 

k— 

oo 

= £>x(fc)-fc(Gi> = A<Gi>. 

k=Q 

In fact, to allow for easy comparison with the data, 
we consider the total conductance per unit area a. By 
making the substitution a = G/A. It follows directly 
that the average value of a satisfies 

{a)=tN A {Gx). (3) 

This is a very intuitive result, since (a) equals the average 
number of dopants in the barrier per unit area multiplied 
by the average conductance per channel. 

To obtain the standard deviation we first observe that 

(Gt)= ((G[ 1) + ... + G[ k) r) = 

= fc(G?) + fc(fc-l)(G 1 ) 2 , 

when the G\ are independently identically distributed. 
Using this fact, we find that 

(G 2 > = f g 2 f(g, A) dg = £ P x (k) J g 2 f k (g) dg = 

k—0 

oo 

= ^P A (fc)-(fc(G 1 )+fc(fc-l)(G 1 ) 2 ) = 

fc=0 

= A(G 2 )+A 2 (G!) 2 . 
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Finally, this yields 

sd(G) = V(G 2 ) (G) 2 = V / A(Gf). 

By making the substitution a = G/A once more, we 
find that the standard deviation sd(<r) (which can be in- 
terpreted as the spread in a) is given by 

sd(a) = ^<G 2 ). (4) 

The most important observation from this equation is 
that sd(cr) is proportional to 1/yA, showing that the 
area-dependence of sd(N)/A leads to a similar behavior 
of spread in a. This also justifies the choice of the fit- 
function in Fig. |21 

APPLICATION TO THE DATA 

Our simple model captures the general features of the 
data and yields reasonable values for the parameters. To 
demonstrate this, a least square fit of the function C j\f~A 
(cf. Eq. 0; C is a fit-parameter) to the standard devia- 
tions has been performed. The result is the dashed line 
in Fig. |2] The fit gives a good description of the data, 
showing that the spread in the data is consistent with 
the prediction of the model. 

From the fit, we make some estimates for the physical 
quantities in the model. The value of the fit-parameter 
C = 1.8 nA V- 1 nm- 1 should be equal to y/tN d (Gf) 
(according to Eq. i jlj l). By looking at the large area val- 
ues in Fig.Hl we find that tN d {Gi) « 0.04 nA V" 1 nm" 2 
(Eq. J2J). Combining these numbers and assuming that 
(Gi) 2 w (G\) we find values for parameters in the model. 
First, 1/tNd w 1/(2500 nm 2 ), which corresponds to an 
average of one dopant per 2500 nm 2 device area (indi- 
cated in Fig.|5J). This number is consistent with e.g. an 
average doping concentration at the interface [lsj around 
Nd = 10 17 cm~ 3 and an effective barrier thickness of 
t = 2.5 nm. Note that because t is the thickness of the 
barrier region where dopants influence the barrier height, 
it is thinner than the total Schottky barrier thickness (a 
few tens of nanometers in this case). Second, we find 
(Gi) ~ 100 nA/V for the average conductance per chan- 
nel. In order to achieve this value, it is necessary to 
have e.g. a small patch with a local barrier height of ap- 
proximately 0.15 eV and an area [lf)j of approximately 
10 nm 2 . These numbers are consistent with the actual 
sample parameters. 

On average, the conductance per unit area increases for 
decreasing island area, as can be seen from Fig. ^ This 
can be explained by the scaling mechanism described in 
Ref. H which predicts that the barrier thickness de- 
creases with device size for devices that are smaller than 
a few times the Debye length. We note that the increas- 
ing value of (a) is not the cause of the increment of sd(cr), 



since similar measurements on low-doped samples do not 
exhibit an increased scatter (not shown). 

As demonstrated by our measurements, the discrete- 
ness of doping is easily observed in highly doped samples. 
In lower doped samples, where the barrier is thicker, the 
effect is expected to be much weaker. The local distor- 
tion of the potential landscape due to the presence of a 
dopant is roughly as large as its effective Bohr-radius, 
which equals about 3 nm in silicon. When the barrier 
thickness is much larger than the dopant's potential well, 
the effective barrier lowering will be negligible. Hence the 
induced local barrier lowering is the most pronounced in 
thin barriers. 

CONCLUSIONS 

In conclusion, we have shown that individual dopant 
atoms dominate the transport characteristics of epitax- 
ial nanometer sized metal semiconductor diodes. Room 
temperature data show increasingly large device-to- 
device conductance fluctuations towards smaller device 
sizes. Measurements at 4.5 K reveal pronounced struc- 
ture in the TV-curves. A statistical analysis based on 
the assumption of randomly positioned individual dopant 
atoms leads to a good description of the experimental 
data. 
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